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ABSTRACT 

We present a new scenario for the cooling of compact stars considering the central source of Cas- 
siopeia A (Cas A). The Gas A observation shows that the central source is a sort of neutron star which 
has high effective temperature, and it is consistent with the well known standard cooling model. The 
observation also gives the mass range of M > 1.5Mq, which is inconsistent with current plausi- 
ble cooling scenario of compact stars. There are some cooled compact stars such as Vela or 3C58, 
which cannot be explained by the standard cooling processes: we invoke some kinds of exotic cool- 
ing processes, where a heavier star cools faster than lighter one. However, the scenario seems to be 
inconsistent with the observation of Cas A. Therefore, we give a new cooling scenario to explain the 
observation of Cas A by constructing models, which include a quark color superconducting phase with 
a large energy gap; this phase appears at ultrahigh density region and reduces neutrino emissivity. 
In our model, a compact star has color superconducting quark core with low neutrino emissivity sur- 
rounded by high emissivity region made by normal quarks. We give cooling curves obtained from the 
evolutionary calculations of compact stars: while heavier stars cool by the standard scenario, lighter 
ones indicate opposite tendency. Furthermore, we suggest that our scenario is consistent with the 
recent observations of the effective temperature during the last ten years. 
Subject headings: dense matter — stars: neutron 



1. INTRODUCTION 

Cooling of compact stars has been disc ussed mainly in 
the context o f neutron stars for decades (jTsuru 13 [1991 
iBeckeii [20091 ). It has been believed that some stars re- 
quire exotic cooling to explain the observed effective tem- 
perature and others can be explained by standard cool- 
ing mechanism, where the central density of the star de- 
termines which cooling process works: an exotic cool- 
ing phase a ppears at higher density above a threshol d 
density fe.g. lYakovlev et all [2001 iGusakov et al.ll2005[ ). 
As a consequence, the heavier star who has high er cen- 
tral d ensity cools faster than lighter one (Lattimcr et al.l 
11991*). However, as described below this scenario be- 
comes inconsistent when we consider the recent observa- 
tion of the effective temperature of Cas A whose mass 
has been found to be unexpectedly large. 

Cas A is the youngest-known supernova remnant in 
the Milky Way and it locates ~ 3.4 kpc distant from the 
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solar system (jReed et al.lll995f ). The supernova explo- 
sion occurred about 330 years ago, but due to absorp- 
tion by interstellar medium, there are no exact historical 
record s except for u nclear detectio n by J. Flams t eed in 
168 (IAshwortbl[l98 0). Recently, IHo fc Reinki ((20091 ) 
and lHeinke fc Hoi (I2010ll have analyzed the X-ray spec- 
tra of Cas A. They give the effective temperature and 
possible regions occupied by mass and radius relations. 
Since Cas A is the isolated remnant, uncertainty of mass- 
radius relation could be large. The lowest mass obtained 
from fitting is about l.SM©. Considering the age 
oi t — 330 yr, Tcs of Cas A must occupy a point of a 
cooling curve due to the standard cooling scenario on 
the (Toff — t) plane. This gives strong constraint on the 
equation of sta te (EoS) an d cooling processes. Further- 
more, .Heinke &: Hd (|2010l ) repo rted the observation o f 
TpfF for Gas A in r ecent 10 years. lYakovlev et all (|2011| ) , 
iPage et al.l (|201lD and iShternin et al.l (|2011[ ) iiisist that 
the rapid decrease in Tee during the years shows that the 
transition to nucleon superfluidity occurs. 

On the other hand, there are some cooled stars whose 
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effective temperature cannot be explained by the stan- 
dard cooling scenario, which does not include the neu- 
trino emission caused by nucleon superfluidity. It needs 
stronger cooling process as in the case of J0205-f6449 
in 3C58 (hereafter, "3C58") or Vela pulsar (B0833-45). 
Also an accreting neutron star SAX J1808 requires strong 
cooling. 3C58 and Vela may be explained by the minimal 
coolin g model which includes nucleon superfluidity (jPagd 
|2009() . However, SA X J1808 needs stron ger cooling than 
the minimal cooling fflein ke et al. 1 120081 ). If we consider 
the strong cooling process according to the conventional 
scenario, their masses should become larger than that of 
Cas A, which may be inconsistent with the mass obser- 
vations of double neutron stars : mass of each ne utron 
star is nearly M ~ 1.4 Mq (e.g. iKaspi et al.ll2006[ ). Iso- 
lated stars should have smaller masses compared with the 
cases of binaries, due to the long standing accretion from 
companions. Although a single EoS must be applied to 
all the compact stars, the existing phase of matter de- 
pends on the density. Therefore, the location of the Cas 
A observation on the (Tog — t) plane becomes very diffi- 
cult to interpret if we believe models with strong cooling 
mechanism to explain all the observations of T^s- 

In this paper, we present models which satisfy both 
cases of Cas A and other cooled stars such as 3C58 
and/or Vela, by considering hybrid stars composed of 
quark matter, hadron matter, and their mixed phase, 
where a characteristic property of color superconduct- 
ing phase is utilized. In addition, we also show cooling 
curves of Cas A for the observations in these ten years 
and indicate for the phase transition to the superfluidity 
to be consistent with the observations. 

2. COOLING CURVE MODELS 

We construct a model which includes both quark- 
hadron mixed phase (MP) and its color superconducting 
(CSC) phase. Considering the first order phase transition 
between hadron and quark phases, it would be plausible 
that both phases coexist and form some kinds of MP. As 
is quite similar t o the "nuclear pasta" phase in the crust 
of a n eutron star (jRavenhall et al.lll983tlHashimoto et al.l 
Il984| ). it has be en shown that MP cou l d form geomet- 
rical s tructures (jMaruvama et al.ll2008D : lYasutake et al.l 
(|2009[) have made EoS of MP under a Wi gner-Seitz 
(hereafter "WS") approximation using MIT Bag model 
for quark phase in finite temperature. In the present 
study, we employ an EoS with the same framework us- 
ing the bag constant B = IQQ MeV fm~^, the coupling 
constant as = 0.2, and the surface tension parameter 
a — AQ MeV fm^^. For hadron phase, we adopt results 
of the Bruekner-Hartr ee-Fock (BHF ) theory inc luding 
hype r ons. A, and E~ (iS chulze et al. 1995; Bald o et al.l 
[19981: IBaldolll999h . Since the BHF results is inappropri- 
ate fo r low-density matt er in the crust, we apply EoS of 
BPS ()Bavm et al.lll97lD for the crust. The EoS gives 
the maximum mass 1.53Mq with the radius 8.6 km, 
and the mass lies within the limits of the observation 
of Cas A. Although our EoS is inconsistent with the 
recent obse rvation of the mass M ^ 2Mq of pulsar 
B1913-hl6 (|Demorest et aLllMol ). we could overcome 
this i ssue by adopting other EoS models feg. lAlford et al.l 
[20051) . 

Using the WS approximation, we obtain a cell radius 
of each phase, and calculate the volume fraction of quark 




Fig. 1. — Volume fractions of quark matter phase having particu- 
lar geometrical structures with a bag cons tant B = 100 MeV fm~^ 
and a coupling co nstant cts = 0.2 IIMaruvama et aLl 120071 : 
lYasutake etaIll2009D . 

matter in MP as seen in Fig. [T] It is difficult to calcu- 
late the neutrino emissivity in MP. Therefore, the vol- 
ume fraction F is multiplied to the original quark neu- 
trino emissivity e^^o dlwamotg ,1980;) ; The total emis- 
sivity by quarks is set to be s^, = Fe„fi. We adopt 
well -known standard cooling process for hadronic mat- 
ter ([Friman k Maxwelll [19791 ) : Modified URCA process 
for the higher density region and Bremsstrahlung process 
for the crust. 

The color superconductivity is the key of the present 
study. There are some sorts of quark pairings such 
as CFL or 2SC according to the degrees of freedom 
of quark flavor and color. It is considered that the 
energy gap of A > 10 MeV is very large compared 
with the temp erature of the center of compact stars , 
Tq ~ keV, (eg. lSchmittll2010l ). Once matter becomes su- 
perconducting, neutrino emissivity must be suppressed 
due to the large energy gap and it could be proportional 
to exp {—A/kBT), where T is the temperature at the rel- 
evant layer and fee is the Boltzmann constant. Therefore, 
in the color superconducting phase with a large energy 
gap (A ^ T), neutrino emiss ivity by quarks is almost 
negligible (jAlford et al. Il2008[ ). We note that we do not 
need to consider which kinds of CSC pairing appears, 
if we adopt a large energy gap. For the neutrino emis- 
sivity . the differenc e between CFL and 2SC is unimpor- 
tant /Berdermannl |2008j). 

We assume that the energy gap in color superconduct- 
ing phase is very large and the phase appears in the den- 
sity region above the threshold density pT. For simplic- 
ity, we use the critical volume fraction Fc instead of px • 
If the matter has a volume fraction F > Fq a,t the density 
in the layer of MP, quarks become color superconducting 
state. 

We select stellar masses of compact stars to be 1.0, 
1.3, and 1.5Mq whose central densities are 1.48, 1.82, 
and 2.67 xlO^ g cm~'^, respectively. We take the crit- 
ical fraction Fc to be 0.1, 0.125, and 0.2 which are 
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Fig. 2. — Cooling curves with color superconducting quark 
phases. Solid, dotted, and dashed lines denote the models with 
the masses of 1.50, 1.32, and 1.03 A/q, respectively. Thick gray 
line on the middle panel denotes a 1.50 A/q model with nucleon 
superfluidity and carbon envelope. Dot-dashed lines with marks 
on the middle panel indicate the model of the mass 1.03 A/q ex- 
cept for the neutrino emissivity in normal quark phase multiplied 
by 1/10 and 1/100 for the lines with marks of circle and triangle, 
respectively. 

appropriate to explain the observations with our sce- 
nario. Results are shown in Fig. [21 with available 
observational values. Most of the observation al data 
are taken from the Table 7.2 in iKaspi et al 
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: ISlane et al.l (I200I) fo r 3C58, IWeisskopf et all (l200l 
for C rab. iPavlov et al I (2001*) for V e la. iBrisken et al.l 
(l2nol for 0656 4-14, FCotthclf eTaP (I2002D for 1706 



44. IKaspi et~all 0UO6) for 1811 - 1925, l^aMereLaLj 



daoda!) for 1823-13, Becker eTall (119961) for 2334-^61. 
Other data are taken from IPaed (I2002D . The data of Cas 
A is attached (IHeinke fc Holl2010ir as the most youngest 
compact star. We find that the cooling curves transit 
from hotter region to cooler region for the parameter 
between 0.1 < Fq < 0.2. As indicated in the middle 
panel of Fig. 121 cooling curves split into two regions 
for Fc = 0.125 and stars with larger masses cool slower 
than those with lighter masses. Since lighter mass stars 
cool faster, they suit for the 3C58 case who does not 
have lower limit of the effective temperature observation. 
However, calculated cooling curves are inconsistent with 
the Vela case which has lower limit and the quark cooling 
is still too strong to explain this case. 

Since the neutrino emissivity of quark phase involves 
large uncertainty, we have calculated the additional cool- 
ing curves for the mass l.O3Af0 in case of Fc = 0.125 
with the neutrino emissivity reduced by a factor of 0.1 
and 0.01. There are some possible factors of this reduc- 
tion for neutrino emissivity accompanying quark /3-decay, 
such as increase of the abundance of strange quarks; de- 
creasing in electron numbers inside MP, l eads to the re- 
duction of the neutrino emissivity (jlwamotg .1980, ). The 



presen ce of 2SC at low density also suppress the emis- 
sivity; iMaruvama et al.l (|2008[ ) discussed that the abun- 
dance of quarks in MP changes and it may cause color 
superconducting phase. We suppose that the reduction 
of the emissivity originates from the above physical pro- 
cesses. If the emissivity of quarks is reduced by these fac- 
tors, the observation of Vela can be explained as shown 
in the middle panel of Fig. [2l 

3. DISCUSSIONS 

We demonstrate the effect of color superconductivity 
in quark-hadron MP on the cooling curve: the larger the 
masses of the compact stars, the slower the speed of the 
cooling. Although the maximum mass of our model is 
at most 1.53 Mq that would be near the lower limit for 
the recent observation of Cas A, our cooling scenario can 
be applied for the cases of M > I.5M0 if more stiff EoS 
is devised. There remain also some problems concerning 
the fundamental physics: uncertain physical properties 
of quark-hadron MP, indefinite tuning of the threshold 
density of CSC, unknown value of the energy gap A. 
Nonetheless, to explain the observed effective tempera- 
ture of Cas A, our model incorporating the color super- 
conductivity with a large energy gap is compatible with 
available observational data. Cooling mechanism associ- 
ated with color superconducting quarks could be plausi- 
ble, because as for the existence of CSC it may be quite 
natural from the rec ent study of phase d iagram between 
quarks and hadrons (jRiister et al.ll2005l ). 

Fundamental physics associated with compact stars 
has still large uncertainty. In the quark phase, the abun- 
dance of each quark is still unknown, and therefore the 
neutrino emissivity is not determined from fractions of 
u, d, and s quar ks and/or c hemical potential of elec- 
trons (Maruyama et al.l 120081) . Although existence of 
CSC phase has been studied well, it is still open to de- 
bate which type of CSC appears, and the quantitative 
values of the critical density and the energy gap should 
be clarified. Since the physical properties of quarks in 
MP are different from those in the uniform phase, there 
would be many factors to change the emissivity. Further 
theoretical study in collaboration with the observations 
is required to constrain the physics of compact stars. 

Considering above uncerta inties, the observat ions of 
the rapid cooling for Cas A (jHeinke fc Holl2010l ) would 
give a clue to constrain some properties of high density 
matter. We have tried to fit the observational data with 
use of a model with nucleon superfluidity. We adopt 
the neutrino emissivity accompanied by the phase tran- 
sition from the normal state to th at of the nucleon su- 
perfluidity (jYakovlev et al.l (|2004l )). where we tune the 
critical temperature of neutron superfluidity. Since 
iHo fc Heinkd (|2009l ) reported that the surface composi- 
tion of Cas A will be carbon and/or helium to reproduce 
the observations of X-ray, we choose the surface compo- 
sition to be carbon and small amount to helium. This 
is because the existence of carbon results in rather high 
effective temperature at the beginning of cooling. To cal- 
culate the cooling curves, we assume the functional form 
of the critical temperature which is phenomenological ex- 
tension of that for the neutron superfluidity as seen 
in the left panel of Fig. [51 This a pproach is similar to the 
method bv lShternin et all (|2011l) except for the profile of 
the critical temperature. Since the cooling curves sensi- 
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tively depend on the critical temperature, which is still 
very uncertain, we draw the cooling curves of upper and 
lower bounds as a band as seen in the right panel of Fig. 
[3l As a result, we suggest that models of cooling curves 
with both CSC and nucleon superfluidity can satisfy the 
observations of the effective temperature of the compact 
star of Cas A during the last ten years. 




Fig. 3. — Left Panel: Density dependence of the critical temper- 
ature of two models of A and IS. Dot-dashed line indicates the crit- 
ical temperature adopted by ^hternin et al. (2011). Right Panel: 
Cooling curves with the nucleon superfluidity with t wo cases of A 
and B which covers the observations of T^g in Cas A I IHeinke fc Hd 
[20T0D . 

There are some experimental projects of hadron col- 
liders with intermediate energy, such as J-PARC or GSI, 
to understand the state of ultrahigh density. They are 
very useful for examining high density character com- 
posed of hyperons, mes ons, transition to no rmal quark 
matter and CSC phase (jAndronic et al.ir2010[ ). However, 
it is still difficult to reproduce the same phase as in the 
core of compact stars by the colliders, where they pro- 
duce still large temperature (T ~ 10 MeV) in the high 
density region of p ^ lO^^gcm"'^ compar ed with the core 
of compact stars (jAndronic et al.l[2009( ). Therefore, it 



is worth while to check theories by both observing as 
many compact stars and comparing theoretical predic- 
tions with observations such as the effective tempera- 
ture. From this viewpoint, further observations of Cas 
A and central sources of other supernova remnants such 
as SN1987A are needed to understand the fundamental 
physics in these extreme conditions. 

Even if quark-hadron MP does not exist, our scenario 
of color superconducting core surrounded by exotic phase 
could be applied to cooled object such as Vela. Consid- 
ering a meson condensed phase, a nucleon superfluidity, 
which reduces the strong neutrino emission by mesons, is 
expected to explain the Vela data consistent with cooling 
curves with some mass ranges. The cooling scenario of 
compact stars affects not only for the cooling of isolated 
stars, but also for the stars in binary systems. There 
are some obse rvational data of X-ray transients (e.g. 
IRutledge et al.l i20Q2) which have gravitational energy 
supply on the surface due to the accretion from compan- 
ions, and X-ray bursts could result due to thermonuclear 
reactions on the surface. These systems are worth while 
to re-examine fro m the point of ex otic cooling as was 
partly done in Yak ovlev et al.l (j2004D for a quiescence pe- 
riod. Since binary systems have larger luminosity and/or 
orbital information such as inclination angle, rotational 
period, and a signal of the gravitational wave, more accu- 
rate mass det ection than isol ated compact stars would be 
possible fe.g. IKramer et all 2006; Wcisb eTg et al.l [20T0I: 
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